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Abstract 

Our  research  program  focuses  on  identifying  reaction  processes  important  for  plasma 
enhancement  of  combustion  systems.  Recommendations  for  temperature  and  energy 
dependencies  of  key  ion-molecule  reaetions  are  made  based  on  the  laboratory  measurements 
described  above.  Thermodynamic  data  such  as  heat  capacity,  entropy,  and  enthalpy  are  not 
available  as  a  function  of  temperature  for  many  of  the  ionic  species  formed  in  these  reactions  and 
are  thus  calculated  using  quantum  meehanical  methods.  The  resulting  kinetic  and 
thermodynamic  data  are  being  used  to  develop  detailed  kinetic  mechanisms  for  incorporation 
into  CFD  models.  In  the  present  contribution,  reaction  rate  measurements  and  thermochemical 
computations  pertaining  to  alkyl  ions,  CnH2n-l-l  -H  (n=l-8)  are  incorporated  into  a  detailed 
combustion  mechanism,  and  calculations  of  the  ignition  delay  time  are  performed.  The  result  of 
the  ignition  delay  time  for  an  isooctane/02/Ar/NO  mixture  to  be  used  in  shock  tube  ignition 
studies  is  presented. 


Introduction 

The  development  of  future  aerospace  vehicles  requires  a  precise  understanding  of  the 
associated  fluid  dynamics  and  chemistry-flowfield  coupling.  Furthermore,  proposed  magneto  gas 
dynamic  (MGD)  applications  require  a  detailed  understanding  of  the  plasma-chemistry-flow 
field  coupling.  While  many  aspects  of  plasma  and  combustion  chemistry  are  known,  the 
coupling  of  these  two  areas  is  not  well  characterized.  This  paper  describes  an  experimental  and 
theoretical  computation  program  aimed  at  developing  a  better  understanding  of  the  important 
and  relevant  ion  chemistry  in  these  plasma  enhanced  combustion  environments. 

Mixing  and  ignition  are  of  fundamental  interest  in  combustion  and  propulsion  research.  It 
is  a  general  assumption  that  for  similar  conditions  the  better  (faster)  the  mixing  and  ignition  the 
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better  the  eombustion  effieieney.  Previous  efforts  have  established  the  effeetiveness  of  utilizing 
high-temperature  plasma  jets  to  eause  ignition  and  rapid  flame  propagation  in  hydroearbon-air 
mixtures.  ’  Renewed  interest  in  the  development  of  high-speed  airbreathing  propulsion 
teehnology  has  prompted  a  number  of  investigators  to  eonsider  adopting  plasmas  as  ignition  aids 
in  seramjet  eombustors  with  positive  results.  However,  eomputational  studies  attempting  to 
assess  the  importanee  of  ion  ehemistry  in  eombustion  have  been  impeded  by  a  laek  of  kinetie 

o 

data  at  relevant  temperatures.  Ionization  and  reeombination  proeesses  at  high  temperatures  have 
been  studied  in  flames  and  shoek  tubes  dating  baek  many  years. However,  the  study  of 
individual  ion-moleeule  reaetions  under  eontrolled  eonditions  has  been  studied  mainly  at  or  near 
room  temperature.  To  date,  numerous  air  plasma  ion  reaetions  with  hydroearbons  at  high 
temperature  have  been  investigated  in  our  laboratory,  and  these  data  are  being  eompiled  for  use 
in  eomputational  modeling  studies  regarding  the  effeets  of  ionization  in  eombustion.  Reaetions 
of  various  atmospherie  plasma  ions  (NO^,  N2^,  Oi  -,  O^,  N^,  HsO^)  with  aliphatie^’^  and  aromatie 
hydroearbons  ’  have  been  studied  as  a  funetion  of  temperature  using  flow-tube  methods. 
Ineorporating  these  ionie  reaetions  into  detailed  neutral  hydroearbon  eombustion  kineties  models 
demonstrates  that  ionization  speeds  the  rate  of  eombustion.  ’  The  detailed  neutral  meehanism 
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used  previously  at  AFRL  is  based  on  the  validated  n-heptane  meehanism  developed  by 

Lindstedt  and  eoworkers  ’  with  the  addition  of  several  reaetions  pertaining  to  the  eombustion 

of  branehed  alkanes  to  model  isooetane  eombustion.  Measurements  of  the  spontaneous  ignition 

of  isooetane  obtained  in  a  rapid  eompression  maehine  at  a  pressure  of  ea.  10  atm  over  the  900- 

95  OK  temperature  range  have  been  well  deseribed  using  the  modified  n-heptane  C/H/N/0 

meehanism.  The  CONP  program  of  the  CHEMKIN  II  software  series  was  used  for  the 

eomputations.  There  is  exeellent  agreement  between  the  ealeulated  ignition  delay  time  and  the 

experimentally  observed  values  as  a  funetion  of  the  eompressed  gas  temperature,  taken  as  the 
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initial  temperature  in  the  eomputations. 

Our  first  studies  foeused  on  evaluating  whether  introdueing  weakly  ionized  air  into  a 
ramjet/seramjet  eombustor  would  signifieantly  improve  hydroearbon  eombustion.  ’  Some 
interesting  eonelusions  of  this  work  are  that  the  residual  ionization  levels  entering  the  eombustor 
from  any  preionization  of  air  exterior  to  the  vehiele  would  not  be  suffieient  to  signifieantly  affeet 
eombustion,  i.e.,  the  ionization  levels  would  be  redueed  to  sub  part-per-million  levels  due  to 
reeombination  before  reaehing  the  eombustor.  These  ealeulations  further  suggest  that  a  plasma 
souree  must  be  present  internal  to  the  eombustor  to  produee  the  ionization  levels  neeessary  to 
signifieantly  reduee  the  ignition  delay  time.  This  eonelusion  is  being  evaluated  in  test  eondueted 
at  the  eontinuous  flow  seramjet  faeility  at  AFRL/PR  at  WPAFB  ’  where  plasma  sourees  are 
being  developed  for  use  inside  the  eombustor  for  ignition. 
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Placing  a  plasma  source  inside  a  supersonic  combustor  results  in  a  new  and  significantly 
more  diffieult  modeling  problem,  because  generating  the  plasma  in  a  hydrocarbon  environment 
greatly  increases  the  number  and  types  of  reactions  that  need  to  be  considered.  Furthermore, 
models,  whieh  hope  to  accurately  predict  plasma  gas  dynamic  effects  in  the  eombustor,  must 
inelude  hydroearbon  ion  chemistry  in  order  to  be  able  to  fully  represent  the  ionization  level  of 
the  gas.  Knowledge  of  such  reactions  is  crucial  for  mechanism  development  and  can  be  a 
daunting  task  for  jet  fuels  consisting  of  mainly  C8-C12  hydrocarbons,  because  the  number  of 
neutral  reactions  that  must  be  ineluded  increases  dramatically  with  increasing  earbon  number  of 
the  fuel.  Luekily,  only  a  relatively  small  number  of  seeondary  hydroearbon  ions  have  been 
encountered  for  the  large  number  of  ion-molecule  reactions  studied.  '  ’  For  example,  for  air 
plasma  ions  reaeting  with  an  alkane  series  given  by  the  formula  CnH2n+2,  ions  with  the  formula 
CmH2m+i^  (m<n)  are  formed  predominately.  Specifically,  for  air  plasma  ions  reacting  with  larger 
alkanes  up  to  octane  (n>5),  C3H7’^,  C4H9^,  CsHu^,  CeHis"^,  C7Hi5^,  ions  are  formed  and 

account  for  a  large  fraction  (ca.  65%)  of  the  reactivity.  Table  1  below  shows  the  ionization 
energy  of  some  simple  hydroearbon  ions  formed  in  these  reactions  compared  with  alkali  ions  and 
eommon  air  plasma  ions.  Values  are  taken  from  the  NIST  database.^®  The  ionization  energies  of 
the  hydroearbon  species  are  significantly  lower  than  the  air  species  and  are  almost  as  low  as  the 
alkali  metals  in  some  eases.  Therefore,  one  may  expect  that  these  ions  are  going  to  be  the 
positive  charge  carriers  in  large  hydrocarbon  combustion  environments.  Consequently,  the  ion- 
molecule  and  eleetron-molecule  reactions  that  determine  the  life  cycle  of  these  ions  will  be  an 
important  consideration  for  any  plasma-based  applieations.  This  statement  follows  from  the  fact 
that  these  ions  may  govern  the  eleetron  lifetimes  due  to  the  large  dissoeiative  reeombination 
rates  observed  for  these  ions.  '  In  most  cases  the  dissociative  recombination  rates  of  these  ions 
is  almost  an  order  of  magnitude  larger  than  that  for  NO+  at  eombustion  temperatures. 

Table  1:  Ionization  energies  (IE)  in  eV  of  some  simple  hydroearbon  speeies.  The  ionization 
energies  of  the  air  plasma  and  alkali  speeies  are  shown  for  comparison. 


Species 

IE  (eV) 

N2 

15.6 

CH4 

12.5 

O2 

I2.I 

NO 

9.26 

C2H3 

8.9 

I-C3H3 

8.7 

i-C3H7 

7.36 

t-C4H9 

6.70 

CsHu 

6.65 

C-C3H3 

6.6 

Na 

5.14 
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Figure  1;  Reaction  path  block  diagram  illustrating  the 
important  ion  processes  important  in  a  plasmaaugmented 
hydrocarbon  combustor 


Because  air  plasma  ions  are  rapidly  converted  to  more  stable  hydrocarbon  ions  in  the 
combustor,  it  is  necessary  to  include  reactions  of  hydrocarbon  ions  with  air  components,  fuel 
components,  combustion  byproducts,  and  electrons.  A  reaction  path  block  diagram  illustrating 
these  points  is  shown  in  Figure  1.  In  particular,  reactions  are  being  studied  with  a  variety  of 
species  including  02,  O,  03,  OH,  H,  alkanes,  alkenes,  alkynes,  and  aromatics  in  our  laboratory 
and  dissociative  recombination  reactions  capable  of  identifying  the  neutral  products  are  being 
conducted  in  collaboration  with  Larsson  and  coworkers  at  the  ion  storage  ring  in  Stockholm. 
Selected  important  and  interesting  aspects  of  these  unpublished,  results  are  outlined  below. 


Secondary  Hydrocarbon  Ion  Reactions 

The  plasma  chemistry  of  ionized  air  reacting  with  hydrocarbons  is  investigated  using  fast 
ion  flow  tube  kinetic  instrumentation.  The  variable  temperature-selected  ion  flow  tube  (VT- 
SIFT)  instrument  can  be  used  to  measure  both  reaction  rate  constants  and  branching  fractions 
as  a  function  of  temperature  from  80  to  550  K.  Measurements  at  higher  temperatures  are 
conducted  using  a  second  fast  flow  tube  instrument,  the  high  temperature  flowing  afterglow 
(HTFA)  instrument.^^’^^  The  HTFA  provides  data  from  300  to  1800  K  enabling  the  study  of  ion- 
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molecule  reactions  at  temperatures  relevant  to  plasma  and  combustion  chemistry.  Both 
instruments  operate  at  ca.  1  Torr  helium  buffer  pressure,  and  the  kinetics  are  observed  over  ca.  1- 
4  ms  reaction  time.  Each  apparatus  has  been  discussed  previously  and  the  following  discussion 
focuses  on  some  recent  results. 

The  oxidation  of  hydrocarbon  ions  in  plasma  augmented  environments  is  of  fundamental 
interest.  Initial  investigations  involve  hydrocarbon  ions  reacting  with  O2,  because  O2  and  N2  are 
the  most  likely  collision  partners  early  in  the  ignition  process.  Since  N2  is  unreactive  and  O2  is 
present  initially  in  such  a  large  concentration  compared  to  the  combustion  byproducts,  even 
limited  reactivity  with  O2  would  be  an  important  process.  For  example,  the  case  of  the  tertiary 
butyl  ion  (t-C4H9’^)  reacting  with  O2, 

t-C4H9+^  O2  ->  C2H50^+  C2H4O,  AH  =  -299  kJ/mol  (la) 

->  C2H3O+  +  C2H6O,  AH  =  -293  kJ/mol  (lb) 

->  HCO^  +  CsHsO,  AH  =  -158kJ/mol(lc) 

->  C2H2O+  +  C2H6O,  AH  =  -67  kJ/mol  (Id) 

has  several  reaction  channels  that  are  exothermic  by  over  150  kJ/mol.  However,  the  rate  constant 
observed  for  this  reaction  is  immeasurably  small  at  300  K  with  an  upper  limit  of 
10'^^  cm^-s'\  which  is  the  current  detection  limit  of  our  apparatus.  CHs"^,  C2H5^,  i-CsH?^, 
sec-C4H9^,  CsHii"^,  and  i-CgHn"^  are  also  unreactive  with  O2  at  300  K,  k  <  5x10'^^  cm^  s'^, 
despite  the  availability  of  reaction  channels  with  exothermic  hies  of  several  hundred  kJ/mol. 

The  next  step  is  to  study  the  reactivity  of  these  ions  with  other  oxidizing  agents  such  as  O 
atoms.  However,  the  quantification  of  O  atoms  as  a  reactant  poses  a  problem  making  the 
experiments  difficult.  Oxidation  reactions  with  ozone  have  been  conducted  in  the  meantime  and 
have  provided  interesting  results.  Oxidation  by  ozone  may  provide  insight  into  the  O  atom 
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reactions  because  the  small  O2-O  bond  energy.  Furthermore,  density  functional  calculations 
indicate  a  mechanism  involving  initial  complex  formation  of  hydrocarbon  cations  with  O3, 
followed  by  O-O2  cleavage,  and  finally  hydrogen  and  alkane  elimination  to  yield  the  highly 
exothermic  products  listed  in  Table  2.  According  to  this  mechanism,  the  O3  and  O  atom  reactions 
may  share  common  intermediates  and  studying  the  O3  oxidation  may  lead  to  some  insight  into 
the  more  difficult  O  atom  studies. 
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O3  Reactant  Concenlralon  (xio”  cm’) 


Figure  2:  Primary  ion  decay  plots  for  CH3^  (squares),  HCO^  (triangles),  C2H5'^  (circles),  s- 
(diamonds),  and  t-C4H9^  (asterisks)  reacting  with  O3.  The  solid  lines  are  nonlinear  least 
squares  fits  to  the  data  performed  to  determine  the  rate  constants.  Experiments  were  performed 
with  a  reactant  gas  composition  of  5%  O3  in  O2  which  is  possible  since  O2  is  unreactive. 

Table  2:  Reaction  rate  constants  for  reactions  of  ozone  at  300  K  measured  with  the  selected 
ion  flow  tube  (SIFT).  The  calculated  collision  rate  constant,  kc,  and  the  measured  rate  constant, 
k,  are  listed  in  italics.  The  reaction  products  and  branching  fractions  are  listed  under  the  rate 
constants.  Energetics  are  taken  from  the  NIST  Chemistry  WebBook.^^  Neutral  products  are  not 
resolved  in  these  experiments.  Brackets  indicate  that  more  than  one  product  channel  could 
contribute  to  the  branching  fraction.  Products  contained  in  parentheses  are  estimated  based  on 
the  mechanism  discussed  in  the  text  involving  O2  loss.  Note  that  the  reactions  are  considerably 
more  exothermic  if  the  reaction  products  in  parenthesis  are  not  dissociated. 


Reaction 

Products 

k;  [kc] 

(10'^  cm^  s’') 

Branching 

Fractions 

AH,  kJ/mol 

HCO+  +  03^ 

H03^  +  CO 

1.2;[1.3] 

>0.98 

-31.5 

CH3+  +  03^ 

HCO+  +  (H2  +  O2) 

1.7;  [1.6] 

0.66 

-410 

H2C0^  +  HO2 

0.16 

-293 

H202^  +  HC0 

0.10 

-311 

02^  +  CH3O 

0.07 

-54 
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H30^  +  CO2 

<0.01 

-1038 

C2H5'^  +  03^ 

HCO^  +(CH4  +O2) 

0.59;[1.3] 

0.70 

-294 

C2H30^  +  (H2+O2) 

0.30 

-392 

S-C3H7^  +  03^ 

C2H30++(CH4+02) 

0.22;[1.2] 

0.64 

-363 

HCO+  +(C2H6+02) 

0.16 

-200 

CH30^+C2H4  +  O2 

\ 

0.14 

-186 

+CH3COOH 

/ 

-670 

C3H30^  +  2H2O 

0.05 

-674 

CsHgO^  +  HO2 

0.01 

-167 

S-C4H9^+  O3 — >• 

HCO++  (C3H8+O2) 

0.08;[1.1] 

Major 

(-0.40) 

-188 

C2H30A(C2H6+02) 

Major 

(-0.35) 

-339 

C2H50^+C2H4+02 

\ 

Minor 

(-0.14) 

-273 

+  CH3COOH 

J 

Minor 

-758 

C2H20^+C2H6+02 

(-0.09) 

-113 

+  HOC2H40H 

I 

Trace 

-423 

CH30^C3H30^ 

C3H50+ 

(<0.02) 

t-C4H9^  +  O3 — >• 

C2H30A(C2H6+02) 

<0.005;[1.1] 

-284 

HC0++(C3H8+02) 

-133 

t-CsHii^  +  03^ 

C2H5CO+  +  (C2H6 

<0.005;[1.1] 

-296 

+02) 

-255 

C2H30A(C3H8+02) 

HCO++(C4Hio+02) 

-112 
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All  of  the  cations  listed  showed  no  reactivity,  k  <  5x10'  cm  s'  ,  with  O2. 

The  following  results  have  been  observed  for  reactions  with  O3.  The  total  reaction  rate 
constants  for  the  alkyl  carbocations  decrease  dramatically  as  the  order  of  the  reactant  carbocation 
increases  as  shown  in  Figure  2  and  listed  in  Table  2.  In  Figure  2,  the  slope  of  the  line  is 
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proportional  to  the  reaction  rate  constant  in  our  experiments.  The  primary  carbocations  CH3+ 
and  C2H5^  react  at  100%  and  46%  of  the  collision  rate  given  by  the  Su-Chesnavich  equation, 
respectively.  The  secondary  s-CsHy^  and  s-C4H9^  carbocations  reacted  with  O3  at  19%  and  7%  of 
the  collision  rate,  respectively.  The  tertiary  carbocations  t-C4H9^  and  t-CsHn^  were  found  to  be 
unreactive  with  O3,  k  <  5x10'  cm  s'  ,  which  is  the  detection  limit  of  our  apparatus  using  this 
ozone  source. 

The  HCO^  cation  is  a  major  product  of  the  primary  and  secondary  carbocation  reactions 
with  O3  listed  in  Table  2.  HCO+  is  known  to  be  unreactive  with  O2,  and  that  finding  was 
confirmed  in  this  study.  The  major  product  of  HCO^  reacting  with  O3  is  H03^.  The  branching 
fraction  reported  in  Table  2  is  >0.98  and  is  intended  to  reflect  the  uncertainty  of  the  source 
conditions.  The  measured  rate  constant  of  1.2  x  10'  cm  s'  is  nearly  equal  to  the  thermal  capture 
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rate  constant  of  1.3  x  10'  cm  s'  . 

The  methyl  cation  was  prepared  by  ionization  of  methyl  bromide  and  was  selectively 
injected  into  the  flow  tube.  However,  despite  very  low  injection  energies,  approximately  12%  of 
the  CH3^  collisionally  dissociated  to  form  CH^  and  CH2^.  While  these  breakup  products  do  not 
interfere  with  the  reaction  rate  constant  determination,  they  do  react  with  O2  to  form  some  of  the 
same  products  observed  in  the  CH3^  +  O3  reaction.^®  Namely,  CH^  and  CH2^  react  with  O2  with 
rate  constants  of  9.7x10'^°  cm^  s''  and  9.1x10'"’  cm^  s'',  respectively,  and  produce  the  following 
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products 

CH^  +  02^  HCO^  +  O  (major),  (2a) 

^O^+HCO  (minor),  (2b) 

^  CO’’'  +  OH  (minor),  (2c) 

and 

CH2^  +  02^  HCO^  +  OH  (major),  (3a) 

^  H2CO+  +  O  (minor).  (3b) 

The  methyl  cation,  on  the  other  hand,  does  not  undergo  two  body  reactions  with  O2  and  only 
exhibits  a  slow  ternary  association  reaction  8.6x10'  cm  s'.  ’ 

Therefore,  the  impurity  CH^  and  CH2^  ions  were  eliminated  before  entering  the  reaction 
zone  by  adding  a  small  amount  (5  SCCM)  of  O2  approximately  20  cm  upstream  of  the  first 
reactant  inlet.  The  and  CO^  ion  resulting  from  reactions  are  converted  to  02^  via  charge 
transfer  with  O2,  and  the  H2CO^  is  converted  to  HCO^  (major)  and  H202^  (minor)  via  reaction 
with  O2.  After  the  addition  of  O2  at  inlet  1,  the  resulting  ion  composition  in  the  flow  tube  at  the 
reactant  inlet  was  observed  to  be  approximately  88%  CH3’’',  8.5%  HCO’’',  2.5%  02^  and  ~1% 
H202^.  The  02^  and  H202^  ions  were  accounted  for  by  subtracting  their  baseline  intensities 
observed  at  zero  O2/O3  reactant  flow  from  the  appropriate  product  ion  intensities.  The  HCO^  ion 
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counts  observed  at  zero  O2/O3  reaetant  flow  were  eorreeted  for  by  subtraeting  the  appropriate 
number  of  eounts  at  eaeh  ozone  eoneentration  based  on  the  reaetivity  of  HCO^  with  O3  listed  in 
Table  2. 

The  methyl  eation  reaets  at  the  eollision  rate  with  O3  produeing  HCO^  as  a  major  produet 
and  H2C0’^,  H202^,  02^,  and  H30^  as  minor  produets.  The  produet  ion  eounts  observed  for 
H2C0'^,  H202^,  and  02^  were  signifieantly  above  the  baekground  eounts  discussed  for  these  ions. 
Therefore  these  ions  are  eonfirmed  produets.  The  H30^  produet,  on  the  other  hand,  either  results 
from  direet  reaetion  of  CH3^  with  O3,  whieh  is  highly  exothermie  but  involves  signifieant 
rearrangement,  or  from  seeondary  reaetions  of  HCO^  with  the  very  low  level  of  H2O  present  in 
the  flow  tube.  Seeondary  ehemistry  ean  be  aeeounted  for  by  extrapolating  the  baseline  eorreeted 
produet  branehing  fraetion  baek  to  zero  reaetant  eoneentration.^^  However,  this  approaeh  did  not 
fully  aeeount  for  the  amount  of  H30^  observed.  Therefore,  H30^  is  listed  as  a  traee  produet  in  the 
reaetion  of  CH3^  +  O3. 

The  ethyl  eation,  prepared  by  ionizing  ethyl  bromide,  reacts  with  O3  to  form  two  primary 
products  HCO^  and  C2H30^.  Because  C2H5^  and  HCO^  are  mass  eoineident  two  separate 
experiments  were  performed  to  determine  the  reaetion  rate  eonstant.  One  experiment  involved 
the  addition  of  CH3CI  to  flow  tube  via  inlet  1.  Methyl  ehloride  reaets  with  HCO^  by  proton 
transfer  to  form  CH4Cr  and  was  not  intended  to  reaet  with  C2H5^.  However,  the  addition  of 
methyl  ehloride  did  eause  a  reduetion  in  the  C2H5'^  signal  and  produeed  several  unidentified  high 
mass  ions  eontaining  ehlorine  so  that  the  amount  of  methyl  ehloride  added  was  sueh  that  only 
about  20%  of  the  original  C2H5^signal  remained.  In  a  seeond  experiment  ^^CCHs^  ion  was 
injeeted  into  the  flow  tube.  Statistieally  equal  amounts  of  '^CHO^  and  ^^CHO^  are  assumed  to  be 
formed,  and  the  intensity  of  '^CHO^  is  used  as  a  measure  of  the  '^CHO^  eontribution  to  the  m/z  = 
30  (sum  of  ^^CCHs^  +  ^^CHO^).  The  reaetion  rate  eonstant  was  determined  by  subtraeting  the 
intensity  monitored  at  m/z  29  (12-CHO^)  from  m/z  =  30  (sum  of  ^^CCHs^  +  '^CHO^)  and 
plotting  the  deeay  of  resulting  ion  signal  as  a  funetion  of  O3  flow.  The  two  different  experiments 
produeed  the  following  results.  Using  the  methyl  ehloride  seavenger  three  determinations 
yielded  6.3x10'^^,  6.5x10'^^  and  5.7x10’'°  em^  s’'.  The  '^CCHs^  experiment  in  a  single 
determination  yielded  5.2  xlO’'°  em^  s’'.  The  reaetion  rate  eonstant  in  the  Table  is  the  average  of 
the  four  measurements. 

The  reaetion  produet  branehing  fraetion  was  estimated  by  injeeting  C  form  of  C2H5 
using  high  O3  flows  where  no  intensity  at  m/z  =  29  remained,  i.e.,  C2H5^  totally  reaeted  and  all  of 
the  HCO"^  was  eonverted  to  H03'^.  Approximately  5%  of  the  C2H5^  dissoeiated  to  C2H3'^  upon 
injeetion  into  the  flow  tube.  The  reaetion  produets  originating  from  C2H3^  were  not  identified  or 
eorreeted  for  and  eonstitute  no  more  than  a  5  pereentage  point  uneertainty,  i.e.,  the  original 
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fraction  of  the  C2H3^  reactant  ion  intensity  assuming  it  yields  only  a  single  product  in  the 
reported  values.  Since  there  are  no  minor  channels,  the  presence  of  C2H3^  has  a  small  effect  on 
the  overall  reported  branching  fractions. 

The  propyl  radical  cation  was  prepared  by  ionizing  i-propyl  bromide  in  the  supersonic 
source.  Approximately  5%  of  the  C3H7^  dissociated  to  C3H5^  upon  injection  into  the  flow  tube. 
Because  C2H30^  is  a  major  product  of  the  sec-propyl  cation  reaction  with  O3  and  is  mass 
coincident  with  the  reaction  rate  constant  was  determined  by  injecting  the  ^^CC2H7^ 

(m/z=44)  form  of  the  reactant  ion.  Given  that  the  expected  C2H30^  product  distribution  is  33% 
C2H30^  (m/z  =43)  and  67%  '^CCH30^  (m/z  =44)  then  the  I(C3H7+)  =  I(m/z=44)  -  2*I(m/z  =43). 
The  reaction  rate  constant  was  determined  using  the  corrected  C3H7^  intensity.  The  product 
branching  fractions  were  determined  from  the  same  13C  experiments  as  follows:  I(C2H30^) 
=3*I('^C2H30+),  I(HCO^)  =  I(^^CHO^)  +  I(^^CHO+),  I(CH30+)  =  I('^CH30^)  +  I('^CH30+)  with 
the  remaining  C3  products  exhibiting  only  a  single  peak.  Three  determinations  were  averaged 
and  yielded  16%  CHO"^,  14%  CH30^,  64%  C2H30'^,  5%  C3H30^,  and  1%  C3H60^.  The  reaction 
products  originating  from  C3H5^  have  not  been  positively  identified  but  can  be  deduced  to  be 
either  CHO^  and/or  €21430"^  based  on  the  observed  product  distributions  for  the  C3H7^  and  s- 
C4H9^  (see  below)  reactant  ions.  Note  that  both  €3147^  and  S-C4H9+  show  a  small  amount  of 
dissociation  to  €3145^  upon  injection  into  the  flow  tube  and  that  only  products  common  to  both 
reaction  systems  can  be  possible  reaction  products  of  C3H5^  and  O3.  The  influence  of  the 
on  the  reported  branching  fractions  has  not  been  corrected  for  but  constitutes  no  more  than  a  5 
percentage  point  uncertainty  (original  fraction  of  the  reactant  ion  intensity  yielding  a  single 
product)  in  the  reported  values  of  the  HCO"^  and  C2H30'^  product  channels. 

Previous  theoretical  studies  at  the  Hartree-Fock  and  MP2  levels  of  theory  on  the  character 
of  the  potential  energy  surface  of  the  propyl  cation  conclude  that  there  exist  two  global 
minima,  that  for  the  secpropyl  cation  and  that  for  comer-protonated  cyclopropane.  These 
calculations  showed  that  a  minimum  does  not  exist  that  corresponds  to  the  npropyl  cation,  which 
appears  to  be  only  a  transition  structure  in  the  interconversion  of  the  sec-propyl  cation  and  the 
corner-protonated  cyclopropane  cation.  According  to  these  calculations,  the  n-propyl  cation  is 
80.5  kJ/mol  higher  and  the  corner-protonated  cyclopropane  cation  is  30  kJ/mol  higher  in  energy 
than  the  sec-propyl  cation  structure.  Furthermore,  experimental  investigations  of  the 
interconversion  of  the  propyl  cation  structures  demonstrate  rapid  rearrangement  to  the  sec-propyl 
cation  form.  A  detailed  study  by  Ausloos  and  coworkers"^^  showed  that  n-propyl  cations 
isomerize  intramolecularly  to  either  isopropyl  ions  or  protonated  cyclopropane  ions  within  10'^^ 
s.  Rearrangement  to  the  isopropyl  ion  is  favored  and  increases  in  importance  with  increasing 
internal  energy  content  of  the  ion.  Mass  spectrometric  investigations  suggest  that  the 
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isomerization  of  the  eorner-protonated  eyelopropane  strueture  to  the  isopropyl  strueture  requires 
10'^  to  10'^  All  of  these  isomerization  times  are  mueh  shorter  than  the  approximately  1  ms 
it  takes  the  ions  to  leave  the  souree  region  and  travel  to  the  reaetant  inlet.  Therefore,  the 
reaetivity  shown  in  Table  2  for  the  propyl  eation  has  been  assigned  to  the  s-CsHy^  isomer. 

To  eonfirm  this  hypothesis,  the  propyl  eation  was  also  prepared  by  ionizing  n-propyl 
bromide  in  the  supersonie  souree.  The  kineties  measurements  were  performed  using  the 
^^CC2H7’^  ion  and  eorreeted  for  produet  as  deseribed  in  the  experiment  using  i-propyl 

bromide  in  the  souree.  The  produet  branehing  distribution  was  determined  from  the  mass 
speetrum  reeorded  with  100  SCCM  of  O3  where  no  remains.  As  above  no  eorreetion  for 

the  produets  originating  from  CsHs^  (~  5%  of  the  original  reaetant  ion  intensity)  was  made.  The 
reaetion  rate  eonstants  and  the  produet  distributions  for  the  C3H7^  ions  generated  from  n-propyl 
bromide  is  found  to  be  nearly  identieal  to  that  produeed  from  i-propyl  bromide.  The  minor 
differenees  observed  in  the  branehing  ratios  ean  be  attributed  to  statistieal  deviations  in  the 
reprodueibility  of  the  measurements  and  small  differenees  in  the  souree  and  injeetion  eonditions. 
This  result  is  eonsistent  with  the  theoretieal  and  experimental  studies  diseussed  above  regarding 
the  stability  and  isomerization  proeesses  of  the  C3H7^  ions,  whieh  indieate  that  only  the  s-  C3H7^ 
isomerie  form  of  the  ion  is  present  in  our  experiments. 

In  pulsed  ion  eyelotron  resonanee  (ICR)  experiments,  Shold  and  Ausloos^^  observed  that 
C4H9^  eations  formed  by  eleetron  impaet  ionization  of  isobutane,  neopentane,  2,2- 
dimethylbutane,  i-butyl  halides,  and  t-butyl  halides  all  have  the  tertiary  strueture.  Fragmentation 
of  n-alkanes,  2-methylbutane,  3-methylpentane,  n-butyl  halides,  and  see-butyl  halides  produee 
both  s-C4H9^  and  t-C4H9^  with  the  s-C4H9'^  surviving  without  rearrangement  for  at  least  0.1  s. 
However,  in  the  ease  of  the  halides,  a  eollision-indueed  isomerization  of  the  s-C4H9^  to  the  t- 
C4H9^  was  found  to  oeeur.  The  ICR  experiments  of  Shold  and  Ausloos  were  eondueted  at 
pressures  of  ea.  10'  Torr  and  observation  times  ranging  from  10'  to  0.5  s.  The  temperature  in 
the  analyzer  eell  was  kept  at  320  K  and  the  eleetron  energy  was  varied  between  10  and  25  eV. 


O3  Reactant  Concentration  (xio'^  cm'^ 
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Figure  3:  Normalized  C4H9^  eounts  (squares)  plotted  vs.  O3  eoneentration.  The  dashed  line 
is  the  result  of  a  single  exponential  fit  of  the  data  to  determine  the  rate  eonstant  aeeording  to 
Equation  (3).  The  result  of  the  single  exponential  fit  is  labeled  ktotai-  The  solid  line  is  the  result  of 
a  double  exponential  fit  of  the  data  whieh  yields  a  fast  eomponent  (ksec)  with  a  rate  eonstant  of  6 
±  4  X  10'  em  s'  and  a  slow  eomponent  (ktert)  with  a  rate  eonstant  <5x10'  cm  s'  . 

Based  on  the  Shold  Ausloos  results,  the  tertbutyl  cation  was  prepared  by  electron  impact 
ionization  of  tert-butyl  chloride  in  the  supersonic  argon  expansion.  The  reported  reaction  rate 
constant  was  measured  by  increasing  the  flow  tube  pressure  to  0.64  Torr  (throttling  the  roots 
pump)  and  provided  a  measured  value  of  <  3.4x10'  cm  s'  which  is  at  the  detection  limit  of  our 
experiments  with  O3  as  a  reactant  and  therefore  represents  an  upper  limit.  No  reactivity  of  t- 
C4H9^  with  O2  was  observed  either,  i.e.,  k  <  3  xlO''^  cm^  s'\ 

The  sec-butyl  cation  was  prepared  by  ionizing  nbutyl  bromide  in  the  supersonic  source. 
Injection  of  C4H9^  into  the  flow  tube  yielded  approximately  93%  €4119^,  5%  and  2% 

C2H5^.  The  electron  impact  ionization  of  n-butyl  bromide  produces  both  s-  C4H9^  and  t-  C4H9^ 
and  collisions  with  halide  molecules  in  the  source  region  can  convert  s-  C4H9^  to  t-  C4H9^.^^  The 
kinetics  plots  with  the  C4H9^  cation  exhibited  curvature  indicating  the  presence  of  more  than  one 
form  of  C4H9^.  Using  a  large  number  of  flow  points  and  fitting  the  data  as  a  double  exponential 
with  a  non-linear  least  squares  analysis  program  provides  an  estimate  of  the  relative  amounts  of 
the  two  types  of  C4H9^  and  a  measure  of  the  rate  constant  for  the  fast  reacting  species.  Figure  3 
shows  the  normalized  counts  of  the  C4H9^  reactant  ion  plotted  as  a  function  of  the  ozone  reactant 
concentration.  The  result  of  a  single  exponential  fit  to  the  data  yielded  the  dashed  line  associated 
with  the  rate  constant  designated  ktotai.  The  single  exponential  fit  reproduces  the  data  very 
poorly,  however,  a  double  exponential  fit  provides  a  much  better  representation  of  the  data.  Free 
fitting  several  sets  of  data  for  the  two  rate  constants  yields  a  fast  component  with  a  rate  constant 
of  8  X  10'  cm  s'  and  a  slow  component  with  a  rate  constant  <  4  x  10'  cm  s'  .  No  reactivity 
of  either  form  of  C4H9^  with  O2  was  observed,  i.e.,  k  <  3x10'^^  cm^  s''.  It  is  assumed  that  the 
slower  reacting  component  is  the  tert-butyl  cation  (ca.  67%)  and  the  faster  reacting  component  is 
sec-butyl  cation  (ca.  33%). 

Additional  insight  into  the  relative  amounts  of  the  different  forms  of  C4H9^  was  obtained 
using  2-methylbutane  as  the  neutral  reactant  in  hydride  transfer  reactions  as  done  by  Shold  and 
Ausloos^^  and  Meot-Ner  and  Field.^^  The  tertiary  and  secondary  forms  of  C4H9^  have  markedly 
different  rates  for  hydride  transfer  from  2-methylbutane  (i-  C5H12).  Based  on  previous  results,®’"'® 
the  hydride-transfer  rate  constants  for  the  following  reactions 
t-  C4H9+  +  (CH3)2CCH2CH3  >  i-C4Hio,  +  t-  C5H11+, 


k298  =  2x\0-\  \  cm3  s-1 


(4a) 
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and 

S-  C4H9^  +  (CH3)2CCH2CH3  >  n-  C4H10  +  t-  C5H11+ 

yt2P5  =  4xlO''°cm^  s'^  (4b) 

differ  by  over  an  order  of  magnitude.  Figure  4  shows  a  plot  of  the  total  C4H9^  reaction  ion 
concentration  as  a  function  of  2-methylbutane  reactant  concentration.  A  biexponential  fit  of  the 
data  plotted  in  Figure  4  yields  two  rate  constants  of  4.9  ±  1.9  x  10‘'°  cm^  s"'  (35%)  and  4.2  ±1.4 
X  10"''  cm^  s"'  (65%).  These  hydride  transfer  rate  constants  are  in  good  agreement  with  those 
reported  previously,  and  this  experiment  confirms  that  approximately  35%  of  the  ions  formed 
under  our  conditions  are  s-  C4H9'^ 


Figure  4:  Normalized  C4H9’^  counts  (squares)  plotted  vs.  2-methylpentane  (i-C5Hi2) 
concentration.  The  solid  line  is  the  result  of  a  double  exponential  fit  of  the  data  which  yields  a 
fast  component  (ksec)  with  a  rate  constant  of  4.9  ±  1.0  x  10‘'°  cm^  s"'  and  a  slow  component  (ktert) 
with  a  rate  constant  4.2  ±  0.3  x  10'^^  cm^  s'\ 

The  product  branching  measurements  are  considerably  more  difficult  than  those  of  the 
reaction  rate  constant  because  most  of  the  C4H9^  is  the  nonreactive  tertiary  form.  Hence  the 
contribution  of  contaminant  ions  produced  on  breakup  during  injection,  C2H5’^  and  are 

considerably  larger  compared  to  the  products  of  the  s-C4H9^  primary  ion.  The  major  products 
observed  are  HCO^  and  C2H30'^,  and  C2H50’^  and  C2H20^  are  minor  products.  As  shown  in 
Table  2,  reaction  of  C2H5^  with  O3  produces  HCO^  and  C2H30'^,  and  C2H5’^  is  unreactive  with 
O2.  Since  C2H20^  and  C2H50^  were  not  observed  in  the  experiments  discussed  above  that 
contain  the  C3H5^  impurity  ion,  these  products  most  likely  originate  from  the  reaction  of  s-C4H9^ 
with  O3.  Given  this  level  of  uncertainty,  the  only  conclusions  that  can  be  drawn  are  that  HCO± 
and  C2H30’^  are  major  products  and  that  C2H20^  and  C2H50^  are  minor  products  produced  in  the 
reaction  of  s-  C4H9^  with  O3. 

The  production  of  C2H50^  in  the  reaction  of  s-C4H9^  with  O3  is  exothermic  and  is  a 
conceivable  product  channel.  The  reaction  for  such  a  process  is  given  below 
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S-  C4H9^  +  O3  >  C3H50^+  CH4  +  O2, 

AH  =  -393  kJ/mol.  (5) 

Since  C3H50^  is  mass  coincident  with  C4H9^,  the  13C  C2H50^  (m/z=58)  form  of  the  reaetant  ion 
was  injected  into  the  flow  tube.  Assuming  a  statistical  C  product  distribution,  approximately 
25%  of  the  C3H50^  produet  ion  produced  should  be  C3H50'^  (m/z=57).  Product  branching 
fraction  measurements  taken  at  high  resolution  resulted  in  barely  deteetable  amounts  of  C3H30^ 
(m/z=55),  '^C3H30’^  (m/z=56),  and  C3H50^  (m/z=57)  are  produced  in  the  reaction.  Trace 
amounts  of  ions  at  m/z=31  and  m/z=32  were  also  observed  and  have  been  assigned  to  CH30^ 
(m/z=31)  and  (m/z=32).  Note  that  the  CH30^,  C3H30'^,  and  CsHsO^  product  ions  could 

also  originate  from  reaction  of  the  C3H5^  breakup  ion  (5%)  reaeting  with  O3.  Therefore,  the 
CH30'^,  C3H30'^,  and  C3H50’^  produet  ions  are  listed  as  traee  produet  ions  in  Table  2. 

The  CsHii^radieal  eation  was  prepared  by  ionizing  2-chloro-2-methylbutane.  The  reported 
reaetion  rate  constant  for  the  CsHn^  ion  formed  reacting  with  O2  is  <  3x10"'^  cm^  s"'  and  with  O3 
is  <  4x10"'^  cm^  s"'.  The  tertiary  form  of  the  CsHn^  eation  is  the  most  stable  form  of  this 
carboeation  and  electron  impact  ionization  of  the  2-chloro-2-methylbutane  ion  precursor  is 
expected  to  produce  this  ion  exclusively.  The  t-  CsHn"^  cation  is  most  stable  cation  studied  in 
this  series,"^^  and,  like  the  exceptionally  stable  t-  C4H9^  cation,  it  is  not  oxidized  by  either  O2  or 
O3  even  though  several  very  exothermic  products  are  available. 

Thermochemical  Calculations 

Building  a  detailed  kinetic  mechanism  requires  an  accurate  knowledge  of  thermoehemieal 
data  or  reliable  predietive  models  for  such  data.  Calcote  and  Gill  have  recently  completed  a 
eomprehensive  compilation  of  the  heat  of  formation,  entropy,  and  speeific  heat  for  ionic  species 
as  a  funetion  of  temperature  to  3000  K.  Although  this  is  one  of  the  most  eomplete  compilations 
to  date,  many  ions  relevant  to  this  work  are  not  included.  Furthermore,  the  thermodynamic 
properties  of  many  of  the  ions  that  are  included  in  the  compilation  are  derived  from  group 
additivity  methods  or  by  comparisons  to  similar  neutral  speeies.  Therefore,  accurate 
thermodynamic  data  for  an  assortment  of  hydroearbon  intermediates  important  to  developing 
meehanisms  for  plasma-enhanced  combustion  are  still  not  available. 

In  order  to  address  this  important  problem,  ab  initio  quantum  mechanieal  calculations^^  are 
being  performed.  These  caleulations  allow  the  estimation  of  the  heats  and  entropies  of  ionic 
species.^^  An  example  will  be  presented  here.  Density  funetional  theory  has  been  applied  to 
elueidate  the  thermoehemistry  of  propargyl  radical  (H2C3H)  and  cyclopropenyl  radical  (C-C3H3), 
for  both  neutrals  and  cations.  The  commonly  used  hybrid  functional,  Becke’s  B3LYP,  was 
utilized  with  a  large  basis  set  denoted  by  6-31 1++G(3df,2p),  whieh  includes  diffuse  functions  on 
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both  the  carbon  and  hydrogen  atoms.  The  neutrals  have  doublet  ground  states,  while  the  cations 
are  both  singlets.  Both  the  neutral  H2C3H  and  cationic  H2C3H^  are  calculated  to  have  C2v 
symmetry.  Neutral  c-  C3H3  is  planar  except  for  one  H  atom  protruding  out  of  plane;  the  point 
group  is  Cl.  The  c-  C3H3^  cation  is  planar  (D3h  symmetry).  The  symmetry  is  important  in 
determining  entropy  and  free  energy  since  the  rotational  symmetry  numbers  are  different,  from  1 
for  Cl  to  6  for  D3h.  The  ionization  energy  of  H2C3H  was  found  to  be  8.72  eV,  which  may  be 
compared  to  the  NIST  database  value  of  8.67  ±  0.02  eV.  The  heat  of  formation  of  H2C3H^  was 
calculated  by  the  atomization  method  and  yielded  AfHion  =1192  kJ/mol,  which  may  be  compared 
with  the  NIST  value  of  1180  kJ/mol.  For  c-  C3H3,  the  calculated  IE  is  6.18  eV,  which  may  be 
compared  to  the  NIST  derived  value  of  6.6  eV.  The  calculated  AfHion  (C-C4F/)  is  1093  kJ/mol, 
which  may  be  compared  with  the  NIST  value  of  1080  kJ/mol.  The  differences  in  the  calculated 
and  experimental  heats  of  formation  are  consistent  with  previous  work,  but  the  discrepancy  in 
the  c-  C3H3  IE  is  still  under  investigation. 

Eor  the  combustion  modeling  described  in  this  paper,  the  thermochemistry  of  the  cations  is 
needed.  A  harmonic  approximation  is  used  to  obtain  vibrational  frequencies  at  the  same  level  of 
theory  as  used  for  the  structures  and  lEs,  described  above.  The  resulting  frequencies  were  scaled 
by  0.989,  recommended  in  the  literature  on  the  basis  of  comparison  with  experiment  for  a  large 
number  of  molecules.  Thermodynamic  quantities  calculated  thusly  compare  well  with  JANAF 
values  (0-2%  up  to  3000  K)  for  molecules  with  known  thermochemistry.  Error  at  high 
temperatures  may  be  due  to  anharmonicity,  which  may  eventually  be  corrected  for  in  a 
semiempirical  manner.  The  calculated  heat  capacities  (Cp)  of  H2C3H^  and  c-C3H3^  differ  by  25% 
at  298  K,  but  at  high  temperatures  become  essentially  indistinguishable,  differing  by  less  than 
0.1%  at  3000  K.  The  enthalpies  likewise  converge  with  temperature,  albeit  at  higher 
temperatures.  The  entropies  for  these  two  isomers  of  C3H3^  diverge  with  temperature  because  of 
the  higher  rotational  symmetry  number  for  c-C3H3^,  and  thus  the  free  energies  also  diverge,  with 
those  for  c-C3H3^  being  higher. 


Figure  5a:  Calculated  heat  capacities  (Cp)  Figure  5b:  Calculated  standard 
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entropies  (So)  of  H2C3H^  and  e-  €3113'^ 
as  a  function  of  temperature  (dots)  and 
fits  to  JANAF  polynomials  (lines). 


of  H2C3H^  and  c-C3H3^  as  a  function 
of  temperature  (dots)  and  fits 
to  JANAF  polynomials  (lines). 


Table  3:  Chemistry  added  to  detailed  C/H/N/0  hydrocarbon/NOx  mechanism  of  Linstedt 
and  Maurice. 
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1.  Measured  at  AFRL  references  7  and  9. 

2.  Estimated  based  on  AFRL  measurements. 

3.  Calculated  based  on  AFRL  measurements  of  similar  systems. 

4.  Reference  6  and  confirmed  in  ALRL  measurements. 

5.  References  27-30.  Rate  coefficients  are  measured,  but  the  products  and  temperature 
dependence  are  estimated. 

6.  Rate  coefficients,  products,  and  temperature  dependence  are  estimated  based  on  other 
measurements. 

7.  Reference  51. 

8.  Included  in  NOx  sub-mechanism  of  Maurice  and  coworkers  and  repeated  here  for  completion. 
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Other  ions  under  eomputational  study  inelude  s-CsHy"^,  t-  C4H9’^,  CeHy^,  i-  CgHn^,  C8H9^, 
CgHii^  and  for  whieh  insuffieient  data  as  a  funetion  of  temperature  are  available  for 

detailed  modeling  applieations.  For  the  smaller  moleeules,  applieation  of  Gaussian-3  (G3)  theory 
is  feasible,  and  should  yield  more  aeeurate  ionization  energies  (IE).  G3  theory  gives  total 
energies  whieh  eorrespond  effeetively  to  ealeulations  at  the  QCISD(T)/GTLarge  level  on 
MP2(full)/6-31G(d)  optimized  geometries,  ineorporating  sealed  HF/6-31G(d)  zero-point 
energies  and  a  “higher-level  eorreetion”  term.  Density  funetional  theory  (DFT)  should  give  more 
aeeurate  heat  eapaeities  and  entropies  beeause  a  mueh  larger  basis  set  is  praetieal. 

The  G3  IE  of  iso-  C3H7  was  found  to  be  7.46  eV.  If  the  G3  method  is  applied  to  the  better 
DFT  geometries  of  the  neutral  and  eation,  this  result  improves  to  7.38  eV.  The  DFT  value  is  7.37 
eV.  These  figures  may  be  eompared  to  the  NIST  database  value  of  7.37  ±  0.02  eV.  These  results 
are  eneouraging  and  suggest  that  lEs  and  thermoehemieal  data  ean  be  sueeessfully  predieted 
within  the  DFT  and  G3  framework  of  model  ehemistries. 

Kinetic  Computations 

In  our  previous  paper,  eomputations  regarding  the  eombustion  of  pure  isooetane  (i-  CgHig), 
pure  ethylbenzene  (CgHio),  and  a  76  vol%  isooetane  /  24  vol%  ethylbenzene  mixture  with  a 
H3O  /NO  /e'  ionization  kineties  seheme  were  diseussed.  Isooetane  and  ethylbenzene  were 
ehosen  beeause  a  mixture  is  a  reasonable  fuel  ignition  surrogate  for  JP-8  fuel.  Some  of  the  key 
reaetions  ineluded  in  modeling  isooetane/ethylbenzene  mixtures  with  a  H30^/N0Ve'  ionization 
seheme  were  also  diseussed.  These  detailed  kinetie  models  —  some  having  hundreds  of 
elementary  reaetion  steps  —  need  validation.  The  eompanion  paper  AIAA-2003-0703  deseribes 
the  potential  to  use  resonant  enhaneed  multi-photon  ionization  (REMPI)  to  ereate  a  large  initial 
eoneentration  of  NO^  and  low  energy  eleetrons  within  a  shoek  tube.  The  goal  of  the  experiment 
is  to  ereate  a  nearly  instantaneous  initial  ionization  eondition  to  study  of  the  effeet  of  ions  and 
low  energy  eleetrons  on  the  ignition  delay  time  of  large  hydroearbon  fuels.  The  experimental 
eonditions  of  the  shoek  tube  measurements  are  an  initial  temperature  of  1000  K  and  pressure  of  1 
atm.  The  gas  eomposition  in  the  measurements  is  80%  Ar,  17.6%  02,  1.4%  isooetane,  and  1% 
NO  and  the  ionization  is  expeeted  to  produee  ea.  1-5  x  10'"^  mole  fraetion  NOVe'  during  the  ea.  6 
ns  laser  pulse  duration. 

Table  3  lists  the  ion-moleeule  proeesses  that  are  eonsidered  for  this  experiment  where  i- 
and  t-  eonnote  the  iso-  and  tert-  isomers,  respeetively,  and  the  absenee  of  a  prefix  indieates  the 
normal  isomer.  The  reaetion  set  that  needs  to  be  eonsidered  in  this  experiment  is  relatively 
straightforward.  NO^  is  unreaetive  with  N2,  Ar,  and  O2  and  reaets  at  the  eollision  rate  with 
isooetane  fuel.^  Under  our  experimental  eonditions,  dissoeiative  reeombination  is  eompetitive 
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with  hydride  transfer.  At  an  intitial  NOVe'  mole  fraction  of  1x10'"^  (8xl0'"^  cm'^)  approximately 
50%  of  the  NO^  ions  undergo  dissociative  recombination  and  50%  undergo  hydride  transfer. 

Isooctane  reacting  with  NO^  produces  a  large  i-  CgHn^  ion  which  can  recombine  with 
electrons  or  undergo  unimolecular  dissociation  at  our  temperatures.  The  primary  product  of  the 
unimolecular  dissociation  is  C4H9^  and  a  tertiary  structure  is  assumed  here  based  on  the 
precursor  ion  and  the  isomerization  processes  discussed  in  the  experimental  section.  The 
unimolecular  dissociation  rate  has  been  studied  over  the  300-500  K  temperature  range  and  only 
at  ca.  0.4  Torr  pressure.  Meot-ner  et  al.  have  studied  the  unimolecular  decomposition  of  Cy 
ions  and  found  that  at  a  few  Torr  pressure  that  the  high-pressure  limit  had  not  been  reached.  A 
theoretical  effort  is  underway  to  try  to  model  this  process.  At  present  a  Lindemann  form  for  the 
pressure  dependence  that  agrees  with  the  available  data  is  assumed  and  is  listed  in  Table  3. 

NO"^  is  rapidly  converted  to  i-CgHiy"^  and  t-C4H9'^  which  then  can  react.  However,  as 
discussed  in  the  experimental  section,  these  ions  are  unreactive  with  Ar,  Ny,  and  O2. 
Furthermore,  since  both  ions  are  tertiary  ions,  their  oxidation  reactions  (if  any)  are  expected  to  be 
too  slow  to  compete  with  dissociative  recombination  and  are  not  included.  The  next  most 
important  process  is  hydride  transfer.  Given  the  focus  on  ignition,  hydride  transfer  reactions  with 
the  primary  fuel  have  been  considered.  Furthermore,  in  the  symmetric  reaction  of  i-  CgHn^  with 
i-  CgHig,  the  reactants  are  the  same  as  the  products,  and  it  is  not  included.  The  t-  C4H9^  ion 
hydride  transfer  reaction  rate  constants  are  slow  and  decrease  with  increasing  temperature.^^’^^  A 
300  K  rate  constant  of  5  x  10"'^  (T/300)  cm^  s'^  is  used  for  the  t-  C4H9^  ion  reacting  with  i- 
CgHig.  The  rate  constant  at  300  K  is  measured,^"^  and  the  temperature  coefficient  is  taken  from 
Meot-Ner  et  al.  based  on  similar  reactions  with  t-C4H9'^.'^®’^^ 

The  ion  cycles  shown  in  Table  3  are  completed  though  neutralization  reactions  with 
thermal  electrons.  The  NO^  dissociative  recombination  reaction  has  been  well  studied  as  a 
function  of  energy  and  vibrational  level.  Many  dissociative  recombination  rates  for  hydrocarbon 
ions  have  been  measured  by  Mitchell  and  coworkers  '  and  the  available  rates  are  included  in 
Table  3.  However,  the  branching  fractions  for  these  reactions  are  unknown.  In  our  previous 
modeling,  equal  branching  between  the  lowest  energy  dissociation  channels,  namely  those 
leading  to  the  production  of-H  atoms  and  -CH3  radicals,  was  assumed.  The  ignition  delay  time 
for  isooctane  is  a  weak  function  of  the  relative  fraction  of  each  species,  -H  or  -CH3,  produced  in 
the  dissociative  recombination.  This  statement  follows  from  our  previous  study  where 
computations  were  performed  involving  selected  atomic  and  radical  species  important  in  chain 
branching  and  propagation  reactions.  The  ignition  delay  time  of  isooctane  was  reduced  as  the 
mole  fraction  of  all  chain  branching  and  propagation  species  was  increased.  In  general,  O,  H, 
OH,  and  CH3  were  found  to  be  equally  effective. 
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Even  though  the  computations  are  relatively  insensitive  to  the  exact  distribution  of 
dissociative  recombination  products,  they  are  sensitive  to  the  total  number  of  dissociation 
products.  Specifically,  if  three  body  dissociation  channels  occur,  especially  those  that  produce 
two  or  more  atomic  or  radical  species,  then  the  presence  of  these  channels  needs  to  be  included. 
Recent  results^^  conducted  in  collaboration  with  Larsson  and  coworkers  at  the  ion  storage  ring  in 
Stockholm  demonstrate  that  for  C2  to  C3  containing  hydrocarbon  ions,  several  dissociation 
channels  are  produced  with  the  majority  being  three-body  channels.  In  light  of  these  new  results, 
the  dissociative  recombination  channels  for  the  i-CgHn^  and  t-C4H9^  ions  have  been  assigned 
50%  involving  -H  +  -CH3  production  and  50%  involving  single  -H  atom  production  as  shown  in 
Table  3. 

The  reaction  sequences  shown  in  Table  3  describe  the  hydrocarbon  ionization  kinetics  as 
well  as  the  main  air  ionization  and  dissociation  steps.  These  reactions  have  been  added  to  the 
detailed  hydrocarbon/NOx  combustion  mechanism  of  Lindstedt  and  coworkers.  '  In  addition, 
a  simple  flame  ionization  mechanism  of  Pedersen  and  Brown^®  has  also  been  added.  Calculated 
thermochemical  values  discussed  in  the  previous  section  are  used  where  available.  Both  the 
detailed  reaction  mechanism  and  themochemical  data  are  available  upon  request. 


Figure  6.  Species  concentration  and  temperature  as  a  function  of  time  for  1.4%  CgHig, 
17.6%  O2,  1%  NO,  80%  Ar  combustion  (II,=1)  at  a  constant  pressure  of  1  atm,  Tinitial=1000  K 
with  (dashed  lines)  and  without  (solid  lines)  ionization. 

Computations  are  performed  for  pure  isooctane/dry  air  combustion  with  a  NOVe' 
ionization  kinetics  scheme.  Presently,  the  CONP  program  of  the  CHEMKIN  II  software  series20 
is  employed,  which  assumes  a  time-dependent  constant  pressure  process.  This  approach, 
although  simple,  does  allow  realistic  studies  of  the  chemistry  and  the  mechanisms  for  flame 
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initiation  to  be  studied  and  eompared  to  shoek-tube  validation  experiment.  The  results  of  several 
CONP  eomputations  are  diseussed  below. 

Results  of  the  eomputations  at  1000  K  for  isooetane  ignition  with  and  without  an  initial 
mole  fraetion  (1x10'"^)  of  NOVe'  are  shown  in  Figure  6.  Additional  initial  eonditions  inelude 
equivalenee  ratio  of  d)  =  1.0  and  initial  pressure  of  P  =  1  atm.  The  ignition  delay  time,  Xig,  is 
defined  as  the  time  at  whieh  the  gas  reaehes  a  temperature  500  K  above  the  initial  value,  Tinitiai, 
and  is  easily  seen  in  the  figure  as  a  nearly  instantaneous  inerease  in  temperature  on  the 
milliseeond  time  seale. 

A  15%  reduction  in  the  flame  induction  time  is  seen  in  Figure  6  for  an  ionization  fraction 
of  10'"^  relative  to  the  case  of  no  ionization  for  this  example.  Larger  reductions  occur  at  lower 
temperatures,  higher  mole  fractions,  and  leaner  fuel  mixtures.  ’  ’  The  ion-molecule  chemistry 

o 

important  for  promoting  ignition  occurs  on  a  very  rapid  timescale  (10'  s)  compared  to  ignition 
(10'  s)  as  shown  in  Figure  7a.  The  i-  CgHis  mole  fraction  is  essentially  constant  on  the  10'  s 
time  scale,  validating  the  earlier  assumption  that  ion  reactions  with  combustion  byproducts  can 
be  neglected. 

In  Figure  7a,  the  NO^  mole  fraction  decreases  rapidly  due  to  dissociative  recombination 
and  hydride  transfer  with  i-  CgHig.  Accordingly,  an  increase  in  the  i-  CgHig^  concentration  is 
observed.  However,  as  i-  CgHig^  is  being  produced,  it  is  rapidly  consumed  through  dissociative 
recombination  and  unimolecular  decomposition  to  produce  t-  C4H9^.  The  t-  C4H9^  ion  is  very 
stable  and  generally  unreactive.  Within  10'^  seconds,  t-  C4H9^  becomes  the  major  positive  charge 
carrier  and  is  eventually  consumed  through  dissociative  recombination.  The  ion  mole  fraction 
decreases  nearly  by  3  orders  of  magnitude  within  10'^  seconds  resulting  in  the  concomitant 
production  of  N,  O,  H,  and  CHg  (Figure  7b)  from  the  dissociation  recombination  reactions 
discussed  above.  No  such  fast  processes  are  present  in  the  thermal  ignition  mechanism  as  can  be 
seen  in  Figure  7c  where  the  O,  H,  CHg,  N  (below  scale)  and  ion/electron  (below  scale)  mole 
fractions  are  orders  of  magnitude  lower  on  the  same  time  scale.  Furthermore,  neither  the  chemi- 
ionization  mechanism  involving  the  reaction  of  CH  with  O  nor  the  thermal  air  ionization 
mechanism  are  active  under  these  conditions  on  these  time  scales. 


Mole  Fraction 
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Figure  7.  Same  conditions  as  in  Figure  6.  (a)  Major  charged  species  profdes  with 
ionization,  (b)  Major  atomic  and  radical  species  profiles  with  ionization,  (c)  Major  atomic  and 
radical  species  profiles  without  ionization. 

Conclusions 

The  fiow  tube  kinetics  measurements  demonstrate  that  ions  found  in  air  plasmas  do  indeed 
react  rapidly  with  aliphatic  fuel  constituents,  leading  to  rapid  decomposition  of  the  alkane 
molecule  and  releasing  reactive  radicals.  These  reactions  are  typically  much  faster  than  the 
reactions  between  neutral  species  involved  in  hydrocarbon  combustion  initiation.  Hydrocarbon 
bond  cleavage,  with  the  associated  radical  formation,  is  the  crucial  first  step  in  hydrocarbon 
combustion.  Combustion  is  generally  initiated  through  either  thermal  decomposition  of  the  fuel 
or  relatively  slow  reaction  between  stable  neutral  species.  The  shock-tube  measurements 
discussed  in  the  companion  paper  AIAA-2003-0703  will  provide  important  validation  data  to 
substantiate  these  findings. 

Models  for  plasma-based  ignition,  piloting,  and  combustion  enhancement  are  useful 
toward  guiding  scramjet  combustor  design  and  development  efforts.  Techniques  for  improving 
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hydrocarbon  fuel  ignition  characteristics  have  been  shown  to  be  enabling  for  viable  hypersonic 
vehicles  powered  by  scramjets.  The  efforts  described  in  this  contribution  will  be  useful  toward 
the  design  of  efficient  combustors  for  hypersonic  propulsion  systems.  Moreover,  the  model 
development  effort  will  facilitate  our  understanding  of  the  detailed  chemistry  involved  in  these 
evolving  technologies.  Finally,  since  the  methods  being  pursued  are  based  on  detailed 
fundamental  kinetics  treatments,  the  models  developed  within  this  program  are  broadly 
applicable  to  a  range  of  diverse  problems  such  as  combustors  for  gas  turbine  engines,  diesel 
engines,  environmentally  clean  combustion,  spark  inhibition,  and  explosion  limits  in  blended 
fuels. 
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